BI-ORTHOGONAL SYSTEMS ON THE UNIT CIRCLE, REGULAR SEMI-CLASSICAL WEIGHTS 

AND THE DISCRETE GARNIER EQUATIONS 



N.S. WITTE 



Abstract. We demonstrate that a system of bi-orthogonal polynomials and their associated functions corre- 
sponding to a regular semi-classical weight on the unit circle constitute a class of general classical solutions to the 
Garnier systems by explicitly constructing its Hamiltonian formulation and showing that it coincides with that of 
a Garnier system. Such systems can also be characterised by recurrence relations of the discrete Painleve type, for 
example in the case with one free deformation variable the system was found to be characterised by a solution to 
the discrete fifth Painleve equation. Here we derive the canonical forms of the multi-variable generalisation of the 
discrete fifth Painleve equation to the Garnier systems, i.e. for arbitrary numbers of deformation variables. 



1. General Structures of Bi-orthogonality 
Consider a formal complex weight w{z) and its Fourier coefficients {iffclfcez defined by 

(1.1) iviz) = y zvkz', IV, = r ^wiQ)Q-\ 

with support on the unit circle where T denotes the unit circle \Q - 1 with C, = e'", 6 e [0, 2tt). The Toeplitz 
determinants constructed from the these Fourier coefficients are related to averages over the unitary group 
U{n) with respect to the Haar measure by the Heine formula tlOl , IfTTl 



1=1 > \ J UO ^0 l<j<k<n 

(1.2) = det[wi-j\j=o....,n-i, n>l, lo[w] = 1. 

Note that w„ + zv-„ in general and consequently the Toeplitz matrix (iX',_y),j=o,...,H-i is not necessarily Her- 
mitian. 

Notwithstanding the fact that many physically interesting quantities are characterised as averages over 
the unitary group U{n) we wish to emphasise another perspective. Intimately connected with such averages 
are systems of bi-orthogonal polynomials on the unit circle which are orthogonal with respect to the weight 
w{z) underlying the U{n) average, in the sense of (|1.4|l and l ll.5|l . Such systems are equivalent to systems of 
bi-orthogonal Laurent polynomials, first introduced by Jones and Thron [24] and studied subsequently by 
125 1 amongst others, as was shown by Hendriksen and van Rossum [17| and Pastro [32J so that all of our 
conclusions apply equally well to these systems. A particular class of weights of great interest is the generic 
or regular semi-classical class which are parameterised by the co-ordinates and residues of singular points 
{zy)^j and {p/l^j respectively (see jl.li for the definition). Some of the relevant properties of this class are 
summarised in a self-contained way in Section 2. The important fact that is relevant here is that systems 
of bi-orthogonal polynomials and their associated functions with such weights have the property that their 
monodromy in the complex spectral variable z is preserved under arbitrary deformations of the singularity 
co-ordinates This fact was first derived in the context of bi-orthogonal polynomial systems on the 

unit circle in [8J, although it was known to be true for systems of orthogonal polynomials on the line due 
to work by Magnus [29 J. This later result was subsequently extended to orthogonal pol5momial systems 
with a certain type of non-generic or degenerate semi-classical weight by Bertola, Eynard and Harnad 
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Il5l . The isomonodromic character of bi-orthogonal Laurent pol5momial systems with such weights has 
also been estabhshed independently by Bertola and Gekhtman [6J. Here we specifically demonstrate that 
bi-orthogonal polynomial systems on the unit circle with regular semi-classical weights constitute a class of 
general classical solutions to the Garnier systems. We carry out this task in Proposition 3.1 by proving that 
the djmamical equations of the system with respect to the singularity co-ordinates are Hamiltonian 

and that this coincides precisely with the Garnier system ^jv = {qj,pj;Kj,Zj}. The reader should note that 
the term classical refers to two distinct notions: semi-classical orthogonal polynomial systems are ones 
which generalise the classical systems such as the Hermite, Laguerre, Jacobi or any member of the Askey 
Table, whilst a classical solution of a Painleve or Garnier equation is a special solution constructible from 
hypergeometric functions and implies a condition on the parameters and boundary/initial conditions. 

The significance of this observation is that one can reverse the usual argument and use the structures 
derived from approximation theory to deduce new results about the integrable system. One particular 
consequence of the identification of the U{n) averages with the Garnier system is their characterisation by 
recurrence relations of the discrete Painleve type. In the case of the simplest U{n) average with one free 
deformation variable the system was found to be characterised by a solution to the discrete fifth Painleve 
equation, see Proposition 14.11 in Section 4. In this section we derive the canonical forms of the higher 
analogues of the discrete fifth Painleve for the Garnier systems, i.e. for the many deformation variable case 
from the approximation theory structures. We give the explicit coupled recurrence relations for the two 
variable Garnier system in Proposition 14.21 and the arbitrary variable recurrence relations in Proposition 
14. 3[ and these constitute our key results. 

We define bi-orthogonal polynomials {(p„{z),(p„{z)}'^^g with respect to the weight w{z) on the unit circle 
by the orthogonality relation 

(1.3) I T— w(C)(/)„,(Q^„(0 = 6„,,„, m,«eZ>o. 

Alternatively one can express this definition in terms of orthogonality with respect to the monomial basis 

(1.4) |^_|_^(g^„(C)C"'. 

(1.5) X^"^^^^'"'^"^^"^\ 

Notwithstanding the notation, (j)„ is not in general equal to the complex conjugate of (pn- The leading 
coefficients of the polynomials are specified by 

(1.6) = z" + A„z"-i + u„z«-2 + . . . + r„, 

K„ 

(1.7) = z" + A„z"-i + p„z"-2 + ... + ?„, 

Kn 

where again A„, p„, f„ are not in general equal to the corresponding complex conjugates of A„, p„, r„ 
respectively. We define a double sequence of r-coefficients by 

(1.8) r„ = , r„ = , « > 1, ro = ro = 1, 

Kn Kn 

which differ slightly from the standard definitions of the reflection or Verblunsky coefficients a„, in that 
«n - —^n+i- The polynomial coefficients introduced above are related by a system of coupled recurrence 
equations, the first two being flOl 

(1.9) kI - kI_^ = (pn{0)MO), A„ - A„-i = r„f„_i. 

We have a extension of the standard results on the existence of orthogonal polynomial systems on the 
unit circle to the bi-orthogonal setting due to Baxter. 

Proposition 1.1 (f3]). The bi-orthogonal system {(p„, (pn]'^^Q exists if and only ifl„ + Ofor all « e N. 
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m < n 


1/k„ 


m = n 





m <n 


1/k„ 


m = n 



(1.10) 



It is a well known result in the theory of Toeplitz determinants 13411 that 

I„+i[w]I„-i[w] 



1 - Vn, n>l. 



{i„[w]r 

Rather than dealing with (p„ we prefer to use the reciprocal polynomial (/)J,(z) defined in terms of the nth 
degree polynomial (pn{z) by 

(1.11) cp'„{z) ~ Z«(^„(l/Z). 

The generating function of the Toeplitz elements, known as the Caratheodory function 

■ z 



(1.12) 



F(z) 



r dC C+: 
JT2raCC-: 



will also feature prominently in our work. The fundamental object identified in the study 1 8J is the 2x2 
matrix 

(1.13) '"??(r?.)' 

^ ^ ' ^ \(p„{z) -e„{z)/w{z))' 

where the associated functions or functions of the second kind are defined by 

(1.14) e„(z)- r ^^w{Q(P„{Q, n>l, eoiz) ^ Mivq + Fiz)], 

Jj IniL, Q-z 



(1.15) 



<,(z) := -z' 



■w{Q(pn{Q, n>\, e*(z) = Ko[iwo - f (z)]. 



Solutions to the orthogonality relations yield the following determinantal and integral representations 
for the polynomials. 



(1.16) 



bn{z) = / det 



Wq 



Wn-l 
I 1 

Wq 



W- 



(1.17) 



Ki^) = y det 
ill 



w. 



n-j 



W„-j-i . 

zJ 

W-„+i z" 



ZV-j+i z' 



Wi 1 



W- 



W-i 

z" 



_ I„[MC)(C-z)] 



— Km 



f„[ K;(Q(l-zC-i)] 
I«N(0] 



The associated ftmctions have representations analogous to J1.16ll.l7|l 



(1.18) 



^e«(z) = :^7^ det 



Wo 



W- 



= z 



W„-Y ... Wn-j-l ... W-i 
y go ■■■ gj ... gn y 
(Wo ... IV-n+l gn 



„ WMQ(l-zC-^)-^] 
I«+i[zi'(C)] 



(1.19) 

where 
(1.20) 
for Izl 7t 1. 



^<,(z) = (-l)""^xy^det 



n-i 



Wl go. 



„,i WMO(C-z)-i] 
I„+i[w(C)] 



ln{z):=2z f ^MO-^" 
Jt 



2raC ^^'1,-z' 



n>0, 



From the general properties of bi-orthogonality we can deduce that the matrix Y„ obeys a difference 
system. 

Proposition 1.2 ( IITOl ). Assuming x„ i= (or equivalently I„ + O)for n ef'i the matrix Y„ satisfies the recurrence 
relation in n 



(1.21) 



Kn+lZ (pn+l{0) 



Yn. 



K„ \(pn+l(0)z K„+i 

Theorem 1.1 ([TO]). The Casoratians of the solutions (p„, (p*„, e„, e* to the above recurrence relations are 



(1.22) 

(1.23) 
(1.24) 



<pn+l{0) 

b„+i{z)e„{z) - e„+i(z)(/)„(z) = 2 z", 

Kn 

(,;,i(z)<,(z) - <,i(z)<|,;(z) = 2^^z«^i, 



(|)„(z)e;(z) + e„(z)(/,;(z) = 2z". 
Under quite general conditions the matrix system Y„ obeys the following spectral differential system. 
Proposition 1.3 ( IIT9ll ,l8l). Assume that the weight satisfies the moment conditions 



(1.25) 



/ 



Then the matrix Yn satisfies the diff'erential relation with respect to the spectral variable z 



f 



(1.26) ^y„:=A„y„ = ^ 



0„(z) + V{z) - ^z0„(z) 



^«+i(0) 



(0) 



0«(z) 



ze;,(z) 



q;(z)-v(z)- 



Kn+l 



0;(z) 



Yn. 



The utility of such a parameterisation of the spectral matrix A„ will be evident when we make the 
specialisation to the regular semi-classical weights. We will refer to 0„, 0,„ 0* , Q* as spectral coefficients. 

The scalar differential equation system corresponding to the above matrix system is specified in the 
following result. 

Proposition 1.4. The components of the matrix Yn satisfy two second-order scalar ordinary differential equations in 
the spectral variable: <pn{z) or e„(z)/ro(z) satisfy 

(1.27) cp'; +p,(p'„+p2(p„=0, 
while <p*„{z) and -e*^{z)/w{z) satisfy 

(1.28) cp^' +p{(P';+pl(P'„=0. 
The coefficients of the scalar second-order differential equations are 

(1.29) 

and 

0„(Q; + y')-0;,(a + V) Kn^l@n 



^ ^ W @'„ IV n 



(1.30) p2{z) = 



W0„ 



Kn W 



[q„ + y _ ^z&n] [q; - y - %f 0;] c/,„,i(o)<^„,i(o)z0„0;, 



W2 



kIW 



and 



and 

(z-ie; + e;:)(Q; - v) - e;(Q;' - v) k„^i e; 



(1.32) p*(z) = 



W&„ Kn ZW 

+ ^ - [^n-v- ^@:,] ^„,i(o)<^„,i(o)z0„0; 



+ 



W2 kIW2 



Proof. The two ordinary differential equations lll.27l l and (|1.28l l follow from the elimination of (/)* and (p„ in 
(|1.261 respectively. □ 

A consequence of the compatibility between the differential relations l|1.26l l and the recurrence relations 
(|1.21l l is the following collection of recurrence relations for the spectral coefficients Q„, Q* , &„, 0* . 

Proposition 1.5 (|j8l). Given the conditions ofProposition U .l\ the spectral coefficients {Q„(z), O* (z), 0„(z), 0* (z))J^^g 
satisfy the following recurrence relations in n 

(1.33) Q„(z) + Q„-i(z) - I + — z) 0„(z) + (« - 1)— = 0, 



n(0) Kn j Z 

n , '^"+1 (\ n, r w , ^'4n+2{0) K„_i(/)„+i(0) (p„+i{0)W{z) 
^•3^ A m\ + ^ U-i(z - □„ z) + — Z0„+i Z Z0„_i z - — — = 0, 

\ ^,,(0) Kn I K„+ic/)„+i(0) K„(p„(0) (p„{0) Z 

o;„.o;,_,„-(i^.^.)e;,,,-„m.o, 

n /'^H+l , ^"+l(0) J.r^> , K„(^„+2(0) K„-i(^„+i(0) K„+i W(z) 

(1-36) — + ^"7^2 (C2„_i(z) - 0„(z)) + ^ — ^z0„^i(z) 2Q„-i(z) + — — = 0, 

\ Kn (p„{0) I K„+i(p„+i{0) K„(pn(0) K„ Z 

(1.37) Q„^i(z) + 0;(z) - f + —z] 0„,i(z) + ^(z0„(z) - 0;(z)) = 0, 

\ 9)1+1 (0) K„+i / Kn 

n ■^s^ n n ^ ''"+2 ^ ^lilll^^^l^^^ , (/)„+i(0)^„+i(0) ^,^_^ ^^^^ n 

(1.38) LJ„(z) - Li„+i(z) + z + 0n+i(z) + 0„(z) z0„(z) = 0, 

Kn+1 \ Kn+l Kn+2 J K„+\K„ Kn Z 

d on\ r\* / \ . r\ I \ h<n+2 </'n+2(0) \ K„+i / \\ ^C^^) n 

(1 .39) Q„^i (z) + Q„ (z) - + — — z 0„^ J (z) - — (z0„ (z) - 0„ (z)) - — — = 0, 

\Kn+l (p„+i(V) J Kn Z 

/I An\ o.^ ^ / ^ , L , <|>«+l(0) (^«+2(0) \ (/)„+l(0)(^„+i(0) K„+i 

(1.40) Q„(z) - Q„^i(z) + — 1 + — z 0„+i(z) + — — - — = 0. 

Kn+1 \ Kn+\ Kn+2 J Kn+\Kn Kn 

The spectral coefficients &n, 0„ are related to their "conjugates" 0* , Q* via a number of functional (or 
recurrence) relations so that one can characterise the system in terms of either set. We will refer to these as 
transition relations. 

Corollary 1.1 (|8J). The spectral coefficients are inter-related through the following equations 
1-41 — z0„ z) 0„-i(z) = -TTTir®" ^) z&„-i{z), 

(p„{0) Kn-l Cj)ni0) K„-i 

(1.42) Q;(z) - ^0;(z) = Q„(z) - ^z0„(z) + n-^^"^ 



^ ^£^,,+2(0) A 

0„+i(z)+ 0„(z) 



(pn+l{0) Kn 



_^ W(z) 



(1.43) Q;(z) + Q„(z) = ^ 

^n+l 

We will require the leading order terms in expansions of (pn{z), (pn{z), e„{z), e;^(z) both inside and outside 
the unit circle. The following Corollary extends the preliminary results reported in jS). 
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Corollary 1.2 ([8]). The bi-orthogonal polynomials (p„{z), <p*„{z) have the following expansions for \z\ < 6 < 1 
1 

(1.44) —(p„{z) = r„ + [r„-i + r„A„-i] z + [r„_2 + r„_iA„-2 + r„fi„-i] + O(z^), 

(1.45) — (/);,(z) = 1 + A„z + p„z^ + v„z^ + 0{z^), 

Kf, 

whilst the associated functions have the expansions 

(1.46) ye„(z) =z"- A„+iz«+i + [A„+iA„+2 - fi„+2] z"^' 

+ [A«+lp«+3 + A„+3p„+2 - V„+3 - A„+lA„+2A„+3]z"+^ + 0(z"+*), 

(1.47) ye-;(z) = f„+iz"^^ + [f„+2 - fn+lXn+l] z"""^ + [fn+3 " f„+2A„+3 - fn+lpn+3 + fn+lK+lK+n] z"^^ + 0(z""'*). 

The large argument expansions \z\ > 6^ > 1 of(p„{z),(()*„{z) are 

(1.48) —(pn{z) = z" + A„z"-i + ^„z"-2 + v„z"-3 + 0(z"-3), 

(1.49) —(pni^) = r„z" + [f„_i + f„A„-i] z"-^ + [f„_2 + f„-iA„-2 + f„fi„-i] z"-^ + 0(z"-3), 
whilst the associated functions have the expansions 

(1.50) ye-„(z) = r„+iz"^ + [r„+2 - ''n+iA„+2]z"^ + [r„+3 - r„+2A„+3 - r„+i/j„+3 + r„+iA„+2A„+3] z"^ + 0(z"^), 

(1.51) y<(z) = l-A„+iZ"l + [A„+2A„+i - ll„+2] Z"2 + [A„+i|U„+3 + A„+3|U„+2 - V„+3 - A„+i A„+2A„+3] z"^ + 0(z"*). 

Proof. Expansions l|1.44b and | |1.49l l are found by differentiating 

(1.52) K„(p„+i{z) = K„+iZ(p„{z) + (/)„+! (0)(/);(z), 
and 

(1.53) (/)„+! (0)(^„+l(z) = K„+iZ"^^(p„+i{z~^) - K„z"(p„{z~^), 

repeatedly, respectively, and setting the argument to zero. Expansion l|1.461 can be found using 

n p:/r^ n ^"^^^ 1 ^n-l{z) - - _ -.^n-liz) - _ _ , (|'«-3(z) 

(1.54) Z = An— + LA„A„-1 - /,i„J + [lU„A„-2 + l-in-lAn - V„ - A„A„-lA„-2\ + n„-4. 

Kn Kii-l Kn-2 Kn-3 

Expansion l|1.50|l is found by making use of 

(1.55) Z(P„{Z) = ^c/)„+i(z) - f ^y(0)</);(z), 

repeatedly. Expansion (|1.47b can be found using the conjugate analogue of the above equation. □ 

2. The Regular Semi-classical Class of Weights 

Of direct relevance to integrable systems is the regular semi-classical class, characterised by a special 
structure of their logarithmic derivatives 

1 d 2V(z) #1 Pi 

w{z) dz W(z) i-^z-Zj 

and its degenerate cases. Here V{z), W(z) are polynomials satisfying the following generic conditions for 
the regular semi-classical class - 

(i) deg (W) > 2, 

(ii) deg (V) < deg (W) = M, 

(iii) the M zeros of W(z), {zi, . . . , zm] are distinct, and 
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(iv) the residues py = lV{zj)IW{zj) ^ Z>o. 
We have the expansion of the denominator in terms of elementary symmetric functions 

M M 

(2.2) W(z) = J](z - Zj) = Y,{-iei[z,, . . .,Zm]z*^-', e, = 1, 



and of the numerator 
(2.3) 



M-l 



M 



lV{z) = (-)'m,[zi, . . . , zm]z**-i-', mo = J^ pj. 

1=0 j=l 

One explicit example, however not the most general form, of such a weight is the generalised Jacobi weight 



M 



(2.4) 



w{z) - ]^(z - ZjY', Pj e C, supp{wdz) = T. 



Lemma 2.1 (fTj/ [281,1^). Let the weight w{z) satisfy the conditions of Proposition 11.31 and w{e^™) = w{l). The 
Caratheodory function (I1.12P satisfies the first order linear ordinary differential equation 

(2.5) W(z)F'(z) = 2y(z)f (z) + U{z), 

where U{z) is a polynomial in z. 

Note that we do not assume one of the singularities is located at the origin and the next result is a variant 
of Proposition 3.1 in |l8|, which did make that assumption. 

Proposition 2.1 (IS)). For regular semi-classical weights (12. 4P , the functions z0„(z), z0* (z), zQ„(z) and zQ* (z) in 
( I1.26P are polynomials of degree deg zQ„(z) = deg zO* (z) = M, deg z0„(z) = deg z0*(z) = M-l, independent 
ofn. 

Because of the assumption zy the following result also differs in detail with the corresponding result 
in (3. 

Proposition 2.2 (f8l). The spectral coefficients have terminating expansions in the interior domain of the unit circle 
about z = with the explicit forms 



(2.6) (-1)"^^--- 



e„(z) = -ucmz ^ + |«eM-i - wiM-i +eM (w + l)A„+i - (n - 1) ^A„-i + j 



+ 0(z). 



(2.7) (-l)^Q„(z) = -neuz ^ + {ucm-i - ^mu-i + em 



{n + l)A„+i -n\A„ + —] 



+ 0(z). 



(2.8) (-l)M^e;(z) = {n + l)eMZ"' + i - (« + l)eM-i + niu-i + 



(«+2) 



rn+2 
f„+l 



An+2 + nA„ 



+ 0(z), 



(2.9) (-1)'^Q;(z) = (« + VfCMZ-' + \ >M-i - (« + l)eM-i + 



1 



(n + l)A„+i + (n + 2) - A„+2)l ) + 0(z), 
Vr„+i /J J 



and in the exterior domain of the unit circle about z = oo with the explicit forms 



(2.10) — e„(z) = {n + l + mo)z^-^ + \ -(n + l)ei-mi + («+2 + mo) 

Km 



■ A„ 



+2 



1 



(2.11) 0„(z) = (1 + -mo)z" ^ + I - ei - -OTi + (n + 1 + mo)A„+i - (n + 2 + wiq) 

(2.12) fcll^e;(z) = -(n + mo)z^-^ + (nej + mj + (fi + 1 + mo)A„+i - (n - 1 + mo) 



+ (« + mo)A4z™-^ + 0(z"-*), 

''n+2" 



H+2 



</)«(0) 



1 



(2.13) Q;(z) = -^moz^ 1 + I ^mi + (n + 1 + mo)A„+i - (m + otq) 



A„ + 



Tn+l 



z*^-2 + 0(z**-3), 
z^-3 + 0(z'^-4), 
z'^-2 + 0(z**-3). 



A„-i + 

rn 



Proof. These expansions follow from the inversion of 111 .261 1, namely the formulae 



(2.14) o'^'-i(0),« 



2^n+i(0)^„^^ = w[-e'„(p„+i +e„+icl)'„j + y[e„+ic/)„ +e„(/)„+i]. 



(2.15) 
(2.16) 

and the expansions of the polynomials and associated functions as given in I|1.44lll.5l1l . □ 

In addition to the coupled equations of Proposition 11.51 and Corollary II. 1[ evaluations of the spectral 
gfficients at the sineular Doints satisfv bilinear relations. 



coefficients at the singular points satisfy bilinear relations. 

Proposition 2.3 (l8j). For ] = !,.. .,M(i.e. Zj + 0) the evaluations of the spectral coefficients satisfy the 
and functional relations 



recurrence 



(2.18) 
(2.19) 

(2.20) 

(2.21) 
(2.22) 



2/ N K„(|)„+2(0) ^ f . . ,,2/ X 

""(^y) = ^ 77^Zye„(z^)0„+l(Zy) + V-'iZj), 

Kn+l(pn+lyJ) 

K„(^„+2(0) 



q;2(z,) = 



Kn+l(pn+l{0) 



^^n-l(Zy) - 



kI (pniO) 



&„{Zj)&:_,{Zj) + V\zj), 



(0) 



4 MO) 



zj&„(zj)@„.^{zj) + V\zj), 



(0)(^„+l(0) 



■Zj@„{Zj)&„(Zj) 



a„{zj) + V{zj) - ^z^0„(zy) 



Q;(z,)-y(z,)-^e;(z,) 



It is these relations that lead directly to one of the pair of coupled discrete PairJeve equations. 
We note the initial members of the sequences of spectral coefficients {0„}J^^q, [®*„}'^^q, {0„))^^g, {□* 
given by 



, are 



(2.23) 

(2.24) 

(2.25) 
(2.26) 



2-t^0o(z) =1V{z)-kIU{z), 
(^i(O) 

2-^— ^z0*(z) = - 2V{z) - kIU{z), 

2(/)i(0)Oo(z) =Kiz[2y(z) - kIU{z)] - k2(/)i(0)U(z), 
2^i(0)zQo(z) = - Kr{2V{z) + Kpiz)] - Kg(^i(0)z!J(z), 



and observe that U{z) defines the initial values for the recurrences of Proposition [LS] 

Hereafter we restore a singularity of the weight at z = (i.e. po 0) in addition to the previous finite 
ones so that the degrees of W, V are augmented by a unit (now M + 1 and M respectively) and the spectral 
coefficients &,,,&*„ and 0„ , O* are polynomials of degree M - 1 and M respectively. In this setting the spectral 
matrix A„ has a partial fraction decomposition 



which define the residue matrices A„ yfor j = 0, . . . , M. We remark that the spectral matrix has singularities 
at z = and z - co irregardless of the locations of singularities of the weight (i.e. zeros of W) due to the 
bi-orthogonality structure. The j-th residue matrix A„ y at the finite singularity Zy is given by 



(2.28) 



^11, j — 



W'(z,) 



Kn+l 



-Ci„{Zj) - V{Zj) + Zj@„{Zj) 



(0) 



<^„+i(0) 



Zj&„{Zj) 



@n{Zj) 



n*„{zj) - v{zj) - ^e;(z,) 



with 2V(zj) = pjW'{zj) for ;' = 1, . . . ,M, while for = the expression is 

(l 



(2.29) A„,o = («-po)(o o"), 
and for the singular point zm+i = oo it is 

(2.30) A, 



In the following section we will draw heavily on partial fraction decompositions of the spectral coefficients 
which imply the summation identities 



M 



(2.31) 
(2.32) 
(2.33) 
(2.34) 



w 



j=0 
M 



/=0 



W'izj) 

W'{Zj) 



0, 



M 



^ Q„(z,) - V{zj) - ^Zj&„{zj) 

j=0 



W'izj) 



^ o^) - v{zj) - ^e;,(z;) 

;=0 



W'(Zj) 



-in + > pi)^ , 

M 

-in + pj), 

M 
j=0 



For simplicity we parameterise the free singularities Zj{t) as arbitrary trajectories with respect to a single 
deformation variable t so that 



M 



(2.35) 



j=o 



where we include / = in the sum for convenience even though zq = 0. 

Proposition 2.4 (|8|). The deformation derivatives of the system ( |L23t with respect to arbitrary deformations of the 
singtdarities Zj are given by 



(2.36) 
where 
(2.37) 



y - R Y - 

dt 



M 

L 

;=0 



Zi- 



An,j 



Z — Z; 



Boo — 



Kn<pn{0) + K„(p„{0) 





K„ 



A particularly important deduction from l|2.36t are the djmamics of the r-coefficients 



M 



(2.38) 



(2.39) 



Q;_j(z,) + V{zj) 



M 
1=0 



W'izj) 



Compatibility of the spectral derivative | |1.26t and the deformation derivative (|2.36b leads to the Schlesinger 
equations for the residue matrices. 

Proposition 2.5 ([8]). The residue matrices satisfy a system of integrable, non-linear partial differential equations, 
the Schlesinger equations, 



(2.40) 



(2.41) 



A; 



\„,j = [Boo, A„,y] + J^^^^^^^ J ^ j = Q,...,M, 



0<k<M 



In the following section we will require a more explicit form for the Schlesinger equations, which were 
first given in [8J. 

Lemma 2.2. The deformation derivatives of the residues of the spectral matrix are given in component form by 

K lie 
(2.42) —Jl—W{zj)A„,^2 = -^Qnizj) 

(/)n+l(0) K„ 



L 

0<k<M 



1 Zj - Zk 

W'{zk) Zj-Zk 



&n{Zk) 



20„(z;) - 2^z^e„(z^) + n^^ 



-0«(2y) 



2Q„(zi) - 2^—Zk@„{zk) + n 

Kn Zk 



and 

n A'X\ W'l. \A </)„+l(0) Kn4>niO) + K„^„(0) 
(2.43) W (Z;)A„;ii = 0„(Z;) 



k*i 

0<k<M 



1 Zj - Zk ( e„{zj) r K„+i 

w'(i) — |0;:(i;) [""(^^) " - — 



-^-^ p„(Zy) + V{Zj) - —Zj@n{Zj) 



O„(z0 - V{zk) Zk@„{zk) + n 

Kji Zk 



o / \ T/^ \ '^"+1 r-. ^ \ , ^^^^^ 
0„(Zy) - V{Zj) Zj&niZj) + n 



for i = 0,...,M. 

Proof. These follow from the Schlesinger equations l|2.40l l and the transition formulae 



(2.44) 

(2.45) 

(/)„+i(0)(^„+i(0) 



W(z,) 

C2„(zy) @„{zj) = 0„{zj) Zj&„{zj) + n , 



Q„(Z;) + V{Z-) - ^zye„(z^) 
Km 



Q„(z,) - y(z,) z/0n(z;) + n 

' -Kn Z, 



'-1 



for i = Q, . . . ,M and where the ratio W{zj)/Zj is interpreted as the limit of zq ^ for / = 0. 

For each singular point Zj the monodromy matrix My is defined by 
(2.46) Yn{zj + Se^™-) = Y„(zy + 6)Mj, 
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and has the classical, triangular structure 

^\-2mp, y i = 0,...,M, 

where Cj is independent of the Zj, and thus t, and n e Z>o but depends on other details of the weight. 



3. The Hamiltonian Formulation and Garnier Systems 
We fix the singularities in canonical position (see 1201 ), i.e. taking them at distinct points 

(3.1) Zo = 0, Zi, . . . , Zn, Zn+i = 1, Zn+2 = oo, 

with exponents po, pi, ... , Pn, Pn+i ~ p, Pn+2 ~ p<x, respectively, so that the number of finite singularities 
is Af + 2. Note that now we have a singularity of the weight at the origin, i.e. po + 0. The denominator 
polynomial for the weight data is given by 

N N+l 

(3.2) W(z) = z(z - 1) [](z - z^) = zY_y-f^^-'z'^N^x-u 

j=l 1=0 

where the elementary symmetric functions of the singularity positions are denoted e;, / = 0, . . . , N + 1 and in 
particular eg = 1/ £n+2 = and e^^+i = H/li ^j- The numerator polynomial is 



(3.3) 2y(z) 1) 7 + ITT + E ^ - E(-) 

;=1 ;=1 '] 1=0 



Z OTjv+i-(, 



where the last relation defines the coefficients mi, I = Q, . . . ,N + 1 and we observe that mo = po + p + L/li pj 
and mjv+i = po^N+i- We can parameterise the upper off-diagonal element of the spectral matrix, i.e. the 
spectral coefficient 0„, which is now of degree N, so that 



(3.4) @niz) = &^\\{z-qr), 0OO = (n + l + mo)— , 

f K„+i 
r=l 

and thus 

r=l ' 

where the poles qr will play the role of canonical co-ordinates and analogue of the sixth Painleve tran- 
scendent. For notational simplicity we will often suppress the n index dependency as we do not discuss 
recurrences in this variable in this section. Furthermore we will assume throughout generic conditions on the 
dependent and independent variables, namely that 

(i) non-coincidence of singular points, Zj + Zk for j i= k and i,k = 0, . . . ,N + 2, so that W{zj) + for 
; = 0,...,JV + 2 

(ii) avoidance by the co-ordinates with each other c\r + qs for r s and r,s = 1,.. .,N implying that 
©', (q,.) + for r = 1, . . . , N and with the fixed singularities qr + Zj for r = 1 . . .,N and = 0, . . . , N + 2 
and consequently W{qr) i= for r = 1, . . . ,N and &„{zj) i= for j = 0, . . . ,N + 2. 
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(3.6) 
(3.7) 

(3.8) 

and 
(3.9) 



Prior to stating our main result we note some summation identities for / = 1, . . . , N 

1 1 W"{zj) 



E 

0<k<N+1 



L 

0<k<N+l 

Pk ^ 

Zj-Zk W'iZj) [W'(Zy)]2 ' 



Zj-Zk 2W'(Zy)' 



L 

0<k<N+l 
k*i 



0«(Zjt) 



^ &niZj) 1 @n{Zj)W'{Zj) 

W{zu)zj-zu W{zj) 2 [W'{zj)Y ' 



1 



l<s<N 



In addition partial fraction expansions imply the following summation identities 

a = 

CT=1, 

&c.[l + LtxZk-L%iqrl a = 2 

= 0,1 
o = 2 , 



(3.10) 

(3.11) 
forr = 1, ...,N, 




N+l 



0„(zy) 




^ " I ©coll + Ljt=l Zfc - Ls*r '?sJ/ £T = 3 



N+l 



(3.12) 



L0n(z,) 
w 



;=0 



W'(Z;) (Zy - (7r)(Z; - <7s) 



0co[l + Zk - Lli qt + qr + qs] - ^,rqt^^. 



a = 0,1,2 

a = 3 . 

(7 = 4 



for r, s = 1, . . . , N, and 



N+l 



w 



&n{Zj) 



W'(Z;) {Zj - qr){zj - qs){zj - qt) 



= -(1 - dt,s)br,sq\ 



®'n{qs) 



iqr - qtWiqs) 



il-5s,r)6tXr 



®'n{qr) 



{qt - qsWiqr) 



{l-5r,t)6s,tq"t 



®'n{qt) 



+ ^r,s^s,tq\ 



,®'n(.qr) 



{qs - qr)W{qt) 
W'{qr) 1 a 



W{qr) 2 &„{qr) qr 



mqr) 

foT r,s,t = 1,. . . ,N with a = 0,1,2,3 whilst for a = 4 an additional term @oo is necessary. 

Proposition 3.1. Assume that generic conditions apply. The dynamics of the bi-orthogonal system is governed by 
the Hamiltonian dynamics of the Gamier system = {qj, pf, Kj, zj} with co-ordinate qr defined above and momenta 

Pr given by 

Q„(t/,) + V{q,) 



(3.14) 

and the Hamiltonian by 



Pr = -- 



W{qr) 



(3.15) 



K 



@n{Zj) ^ W{qr) 1 



' ~ W'Cz,~) 2_i el 



y^'{^i)j^®'n{qr)Zj-qr 



= Ai^i{qr), r = l,...,N, 
(2V{qr) n 1 \ n(l + mo) 



\W{qr) qr Zj-qr) qr{qr-l)_ 
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The indicia! exponents are given by 6j - -pjfor j - 1, . . . ,N + 1, 6q - n - po, and Uoo = -n, doc - n + 1 + E^o^ Pj 
with the constant k = -n{l + mo). The latter relation is the necessary condition for a classical solution to the Gamier 
system, and for the "seed" solution ofn = it vanishes. 

Proof. We seek to compare our system with the second order ODE for the Garnier system Eqs. (4.1.1) and 
(4.1.4) in EOl, see also S, 1221, 



Using Prop. I1.4l we find that 



Z;(zj - l)Kj yi q,{q, - 
.^^ z{z - l){z - Zj) ^ Lz{z-l){z -qr) 



= 0. 



(3.17) 



po + l-n p + 1 



1 ^ y. + 1 y. 1 



From the residues of pi at z = Zj for - 0, 1, . . . , N + 1 we can read off the indicial exponents. We also note 
that 

(3.18) P2 = [-n{l + mo) + 0{z-')] , 

as z —> oo. 

An alternative form of p2 to | |1.30l l can be given in terms of the residue matrices 



(3.19) 



N+l N 



An 



.J / J ('2 Z Mz £7 ) ^ ^ ' 

j=0 r=l ^ 1'^ 0<j<k<N+-l 



TrAnjTrAnk - TrA„jAnk - Anj,ii - A„,t,ii 



(Z - Z;)(Z - Zfc) 



From p,- = Res2=q,p2(z), the above relation and l|2.28|l we find (13.141 1. We need to invert this relationship and 
express Q„ + V in terms of the canonical co-ordinate and momenta. Because this variable is a polynomial 
of degree N + 1 in the spectral variable z we can use the Lagrange interpolation formula at the nodes 
z = 0,qi, . . . , qN, 1 (which are assumed to be distinct) 



(3.20) n„ + V^J^ [Q„{q,.) + V{q,.)] 



r=l 



z(z - 1) YlsA^ - Rs) 
qAqr - 1) XisA'ir - qs) 

(1 - z) WAz - qs) z WAz - fl,) 

+ [Q„(o) + vm - — /V , ^ + [a(i) + v{i)] ■ '^'^^ 



Ylsi-qs) 



ns(i-^s) 



The coefficients of the two last terms are known as Qn(0) + V{Qi) = {—)^{n - po)eN+i from l|2.7|l and utilising 
the coefficient of z^^'^ in (|2lTt we deduce 



(3.21) 



Q„(l) + y(l) l + mo 



0«(1) 

Consequently we conclude that 



0c. 



+ (-)~(n-po) 



^N+l 



E 



PrW{q,) 



(3.22) 



0„(z) + Viz) ■ 



z@„{z) = @„{z) 



0n(O) ^ qriqr - l)&'n{qr) 

Z PrW{qr) 



N 



'@J ^ ^ ^" ^"^©..(O) ^^z-qrqr@'„{qr) 



E 



We also require a similar representation of 2V(z) and proceeding in the same manner we find 



(3.23) 



2y(z) = e„(z) 



W'(0), W'(l) 



0„(O)' 



0«(i) 



^ 1 

<-»E74 



2V{q,.) 



^z-q,. qriqr - l)@'n{qr) 



From this, under the limit z ^ 00, we have the summation 



(3.24) 



E 



2V{q,) 



_ mo_ W'jO) _ W'jl) 
^ qriqr - l)@'„iq,) " ©o. ^ ^° ©„(0) ^ 0„(1) ' 
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and for z = Zy we have 



(3.25) 



2V{q,.) 



W'iO) 



^ q,.{qr - l)&'niqr) Zj - q,- ^" e„(0) ^' ^ 0„(1) ' e„(zy) 
A number of other sums for j = 1, . . . , N which will be required subsequently are 

2V{q,) mo 2V'{zj) 2V{zj)&„{zj) 



(3.26) 

(3.27) 
and, for r ■ 



T,^{Zj-qrf@'„{qr) 0. 



IViqr) 



N 

y 

(Z/ - qr)qr&n{qr) 



&n{Zj) 

IV'iO) 



®l{Zj) 

IVizj) 



Zj@M Zj@n{Zj)' 



1, . . . ,N the sum 

N 



(3.28) qr{qr-l)Y^ 



2V{qs) 



qsiqs - i)@'„{qs) qr - qs 



W'(0) 



w'(i) 



@n{or' ' ^0«(i) 

The analogous representation for W(z) takes the form 



V'{qr) 1 0;'(<?.) Iqr - 1 



(3.29) 



W(z) = z(z - l)0„(z) 



11 
0l "^I^z^ 



2V(r?r) 
W(^,) 



1) 



0°° ^ " ^Mr - l)0n('?» ) 



which enables one to infer the summations 



(3.30) 
and for r : 
(3.31) 



E 

/-=! 



W{qr) 



1 



{Zj-qr)qr{qr-i)&'n{qr) ©o. 



1,. . .,Nthe sum 

N 

L 



W{qs) 



^ qsiqs - l)0,',('?s) <?r 



1 ^ W(r?,) 

0OO (Jr((?r - l)&'n{qr) 



i = i,...,N, 



W'iqr) 1 &:;{qr) 2qr - 1 



W{q,) 2 @'Mr) qriqr - 1) 



We proceed in two steps - the first being to compute the derivatives of the canonical variables using the 
theory from fS) and the second to compute the Hamiltonian itself. Then all we require is a verification of 
the Hamilton equations of motion. 

Our first task is the computation of the deformation derivative of qr 



(3.32) 



k=l 



. dq Ai2 
= -Res2=q,.-— , 
dzi: An 



and where we use the partial fraction expansion \2.2n\ for the 1, 2 element and the Schlesinger equation 
(|2.42)l for An]xi- Now we employ the expressions for 0„(z) and Qn(z) - ^^^z0„(z) atz - Zj,Zk in terms of 
the canonical variables as given by | |3.4| | and | |3.22| | along with | |3.23| | in this formula. After interchanging 
the summation order each term contains a factor which is a sum over the singularity locations and we can 
evaluate these using the a = Q cases of (|3.10t , l|3.11t and l l3.12t . Simplifying we find for r, j = 1, . . . ,N 



(3.33) 



(zj - q,-) 



dZi 



0„(Zy)W(^,) 
@',MrW{Zj) [ 



2pr + 



2V{qr) 
Wiqr) 



n 

qr 



Zj - qr 



We now proceed to the computation of p, which is a more laborious task. We start with the partial 
fraction expansion for A„ ii(<jr) given by ^2.27} , differentiate with respect to t and employ the 1, 1 component 
of the Schlesinger equation (|2.431 for Anj^n- Into this expression we must substitute the expressions for 
0„(z) and Q„(z) ± V{z) — ^^^z0„(z) at z = Zy,Zs; in terms of the canonical variables as given by l|3.4t and 
l|3.221 along with (|3.23l l. Again we interchange the summation order which yields in each term a factor of 
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a singularity sum. To evaluate these sums we employ all of the formulae given in | |3.10|| , l|3.111 , l|3.12|l and 
(|3.13t and assemble all the terms together. We also require the formula for q,- implied by I I3.33II . Considerable 
cancellation occurs at this stage but we are not yet finished. To arrive at the final result we must employ the 
transcendent sums | |3.24| | and I I3.28I I, and the result for r, / = 1, . . . , N is 



(3.34) (z,-^,,)— 



W'{z,) d W{qr) 



@„{zj) dzj ©;(<?,) f'- \ W{q,.) 2 ©;,(<?,.) ; " W{qr) \ V{q,) 2 ©;,(<?,•) 
TTTTTT - o Z^TTTT " " " Z T TTTTTT " o Z^TTTT + Z T \ + n{l + WZq)- 



W'iq,) l&'^{qr)\ , 2Viq,) (V'{q,) 1 0^'(^,) 



qr \ W{qr) 2 @'„{qr) qr j Zj - qr \ V^{qr) 2 @'„{qr) Zj - qr j qMr - 1)(Z; - qr) 



mqs) 1 

®'n{qs)qs-qr 



^2 



2V{qs) ps ps , qs- qr 



+ Ps - n— V n{l + mo) 



W(qs) qs Zj- qs qs(qs - 1)(Z; - qs) 



Our starting point for the computation of the Hamiltonian is the formula 

N ^ 1 

(3.35) Kj = -A„j,iiy y [TrA„;TrA„fc-TrA„^v4„t-A„y,ii-A„fc,iil, i = l,...,N. 

r=l 0<k<N+l ' 

We substitute (|2.28l l and (|2.29l l for the residues of the spectral matrix in the above formula in conjunction 
with the transition relations | |2.44| | and ll2.45l l. Into the resulting expression we must further substitute 
the representations of the spectral coefficients | |3.4| |, ( |3.22| | and | |3.23l l. Some care needs to be exercised to 
ensure that the j,k = contributions are correctly accounted for and this leads to a separation between 
these terms and the remaining generic ones. In the group of generic terms there is a summation over Zj; for 
< k < N + 1 and in these terms we can employ the evaluations given in I l3.6b , | |3.7|I , | |3.8|I , (|3.101 , B.lHl and 
(|3.121 . After reinstating these terms in the total expression considerable cancellation is effected. Still further 
simplification is possible by using the transcendent sums I |3.24| |, | |3.26l l, i3.27\i , | |3.25l l and lastly | |3.30l l. The 
final result yields | |3.15l l which is precisely Eq.(4.3.7) in |20| after making all the notational correspondences. 
The verification of | |3.33| | and | |3.34| | using this Hamiltonian is a straight forward calculation. □ 

Remark 3.1. The Riemann-Papperitz symbol (see subsection 1.4 of [20J, "Fuchsian Equations" for the defi- 
nition) for our system is 

(Zo = Zi • ■ • Zn Zn+1 = 1 Ziv+2 = oo q^ ■■■ 
0---0 -n ■■■ 

n-po -pi ■■■ -pN -p 1 + Lf=o Pi 2 ■ • • 

Remark 3.2. One could formulate the Hamiltonian system using the "conjugate" set of variables ©* , Q* 
instead, and in fact all of the dynamical description, however we do not carry out this task as no new 
understanding would result from it. 

Remark 3.3. This Hamiltonian system is not polynomial in the canonical variables qr and the dynamical 
equations for qr do not possess the Painleve property in zy, however using the well-known canonical 
transformation to the new Hamiltonian system Hn - {Qj, Pf, Hj, ti) 123/ Il20l 



(3.37) tj = 



2y 



z,-l' 



Z;0„(Z;) 



I' 

N 



(3.39) Pj = -{zj-l)Y^ 



©coW(^,) 



qAqr - I) {qr - Zj)e'„{qry 
both these deficiencies can be removed. 

Remark 3.4. The definitions of the co-ordinates (|3.41 and l|3.141 , l|3.151 bear some resemblance to those in 
Syklanin's separation of variables procedure I15J , LI 6.1 although we cannot offer any explanation of this here. 
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4. The discrete Garnier Equations for the L(1 ; 2) Garnier Systems 

In this section we derive recurrences in n or equivalently when ±1 and 0oo ^ 6oa ± 1 for the 

appropriate variables, thereby characterising the system in this way. We will find the discrete fifth Painleve 
equation as the simplest case corresponding to one free deformation variable and the higher analogues of 
this recurrence system for the multi-variable generalisations. These are the discrete Garnier equations. 

The system with M = 3, N = 1 singularities at the standard positions 



(4.1) 



Oil CO 

n- po -pt -pi n + po + pt + pi 

corresponds to sixth Painleve isomonodromic system and was treated extensively in |7l, especially with 
regard to the forms of the discrete Painleve equations. The weight data is 

(4.2) W(z) = z(z - t){z - 1) = z^ - eiz^ + £'2Z - 63, 2V{z) Y = nioz^ - niiz + nii. 

Z — Z; 
/=0,f,l ■' 

The spectral coefficients can be parameterised in the following way 

(4.3) — e„(z) = 9„ + (n + 1 + OTo)z, 

Kyi 

111 

(4.4) Q„(z) = -ne2 + -m2 - (w,, - -mi)z + (1 + -nio)2^, 

whilst the sub-leading coefficients can be related to the polynomial coefficients themselves by 

(4.5) 9„ = - — {n-po)t, 

^n+l 



(4.6) a)„ = 1 + po + Pt + (1 + po + pi)i + {n + 2 + mo) 



, rn+2 

Art+2 

rn+1 



(« + l + mo)A„+i. 



Proposition 4.1 ((7j). The n-recurrencefor the bi-orthogonal polynomial system with theM = 3 regular semi-classical 
weight is governed by the system of coupled first order difference equations 

^ [co„+n-t- po(f + 1)- {pt + pi)t][cv„ + n-t- po(f + 1) - p, - pit] 
Jnjn+i - ^^^^ p^(^f + 1) _ _ pj][a,„ + nt-l- po{t + l)-pt- pit]' 

and 

1-t 1-t 
(4.8) cOn + co„-i + {In -l)t-2- 2po(f + 1) - 2pf - pi{t + l) = {n- po)- + (« + 1 + mo) 



fn-1 ' 'tfn-l 

The transformations relating these variables to the bi-orthogonal system are given by 
(4 9) ff.-M) 

Proof. We refer the reader to the proof of the next case. Proposition 14.21 as the methods employed in both 
cases are the same. □ 

The above coupled recurrence system is equivalent to the canonical "discrete fifth Painleve equation" 
1131 . Il30l with the mapping 

(4.10) fi-^l/f, {l-t)ah,-nt + l + po{t + l) + pt + pi, 
and the identification 

(4.11) ao = pt, ai=n-po, a2 = -n-pt-pi, as = pi, = n + 1 + po + pt + pi- 

In the case M = 4, N = 2 we have the two-variable generalisation of the sixth Painleve equation or the 
L(l^; 2) two-variable Garnier system. A standard placement of the singularities, without loss of generality, 
would be 



(4.12) 



s t 1 00 

n- po -ps -pt -pi n + po + ps + pt + pi 
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We write for notational convenience the weight data in the following way 

(4.13) W(z) = z(2 - l)(z - s)(z -0 = 2*- eiz^ + e^z^ - ej,z + e^, 

(4.14) 0,7/.X_T.7 V Pi 



2V{z) = W 



Z - Z; 
]=0,s,t,l ' 



otqz^ - niiz^ + W222 - nij,. 



The Toeplitz elements satisfy the third order linear difference equation 

(4.15) (y - po)stWj - [(; - 1 - po)(s + f + st) - pist - pst - pts] Wj-i 

+ [(;■ - 2 - po)(l + s + f) - pi(s + - Ps(f + 1) - Pf(s + 1)] Wj-2 

- (/ - 3 - po - Pi - Ps - pO^^i-s = 0, / e Z, 

and we take lO-i, Wq, W\ as defining a solution of the system. Also, according to ll2.10l l- (|2.7b , we can parame- 
terise the spectral coefficients as 



(4.16) 



(4.17) 



^^0„(z) = {n- po)e3—^ + 9„z + {n + l + mo)z^, 

11 11 

0„(z) = nea - -mj + (w,, - -m2)z + (cD„ + -JfzOz^ + (1 + -mo)z^ 



in terms of 9„, co„ and (D„. 

We can relate the new variables introduced above to the coefficients of the polynomials through the 
explicit forms of the spectral coefficients | |2.10|| - (|2.7|I in the following way 



(4.18) 
(4.19) 
(4.20) 



VVn+l 



dn = -{n + l)ei - OTi + (n + 2 + mq) 

(On = -ne2 + m2 + (n - po)e3 
cD„ = -ei - )«i + (n + 2 + mo) 



-(n + mo)r„+ir„ - 2A„+i, 



Irn+i 
r„+2 
rn+i 



rn+2rn+l 



A 



n+l- 



These formulae only serve to allow the recovery of the original variables and do not feature in the recurrence 
relations. 

At this point it is possible to use the foregoing results to derive a system of recurrence relations which is 
the analogue of the fifth discrete Painleve equation for the two-variable Gamier system. 

Proposition 4.2. The following system of coupled first order recurrence relations in nfor the variables {f„ , g„, cOn, l^g 
completely characterises the bi-orthogonal polynomial system 



(4.21) s/„/„+i = 

(4.22) tgngn^i = 



[a)n + SCO,, + (1 + mo)s^ + {n- po)t + ps{s - t){l - s)] [w,, + scD„ + (1 + mo)s^ + (« - po)i] 
[co„ + cD„ + 1 + mo + {n - po)st + pi(l - s){t - 1)] [a}„ + ©„ + 1 + mo + (n - po)st] 

[aj„ + t(D„ + (1 + mo)P + {n- po)s + p,(i - s)(l - f)] [w„ + to„ + (1 + mo)t^ + {n- po)s] 
[co„ + a)„ + 1 + mo + (m - po)st + pi(l - s){t - 1)] [a)„ + cD„ + 1 + mo + (n - po)st] 



(4.23) co„ + co„-i =m2-{n- l)(s + t + st) 



^ S^-t^ + {l-S^)tg„-{l-t^)sfi, t-S + {l-t)f,-{l-s)gn 

(" " Po) — I — T-t; — rr? M — — + {n + l + mo)st 



t-S + {l-t)fn-{l-s)gn 



(4.24) a)„ + cD„_i = -mi + (« - 1)(1 + s + f) 



^, t-s + (l-f)s/„-(l-s)^g„ ^ ^ 

+ (" " Po)-;^ — 7-t; — TT? — — ir:^ + (n + 1 + mo) 



t-S + {l-t)f„-{l-s)gn 



t - S + {1 - t)sfn - (1 - S)tg„' 

S^ - t^ + {1 - S^)tg„ - {1 - t^)sf„ 
t-s + {l-t)sf„-{l-s)tg„ ' 



where 
(4.25) 



e„(s) 



e„(i)' 



@nit) 

©«(!)■ 
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The recurrence relations are subject to the initial conditions 

(1 + mo)sro_i - [pos(f + 1) + pist + pst + p(s] wq- {1- po)stwi 



(4.26) /o 

(4.27) ^0 



(1 + w2o)ro-i - [po(s + + Pisi + Pst + Pfs] rfo - (1 - Po)sfJi'i 
(1 + mo)tW-i - [pot{s + 1) + pist + Pst + pfs] Wq- {1- po)stoi 
(1 + mo)zy_i - [po(s + + Pist + pst + P(s] toq - (1 - po)stoi 



(4.28) ojo = {1- po)st— + post— + po(s + i + sf) + pist + pst + pts, 

Wq W-i 

(4.29) cDo = -(1 + mo)— - (1 - po)st— - [po(s + t + st) + pist + p^f + pts] — . 

Proof. The first two relations ||4.21I4.22| | follow from the evaluation of (|2.181 at = s,t,l and taking ratios. 
From the definitions l|4.25l l and ||4.16|I we note that 

(4.30) /„ - '""1 



(4.31) g„ 



(m - po)stj^ + d„ + n + l + mo ' 
(" - Po)^T^ + ^" + (" + 1 + mo)t 



(n - po)st-^ + d„ + n + l + mo 

Solving for d„ and r„/r„+i in terms of /,„ g„ we find that 

- F + (1 - s2)te, - (1 - t^)sf,, 
t-s + (l-Os/„-(l-s)tg„ 

and 

r„ n + l + nio t - s + {1 - t)f„ - (1 - s)g„ 



(4.33) - 

r„+i n-po t - s + {1 - t)sfn - {1 - s)tg„' 

The second pair of relations l|4.23l4.24|l follow from the evaluation of the recurrence lll.33l l at z = s, t, then 
solving for their left-hand sides and utilising the preceding expressions. □ 

We now give the result for an arbitrary number of variables in a canonical placement expressed by the 
abbreviated symbol 

f ti ■■■ In 1 oo 

(4.34) \ ^N+i 

[n-po -pi ■■■ -pN -p n + Zj=Q pj 

with M = N + 2. To begin with we denote the Z-th elementary symmetric function in the variables ti,...,t„ 
by e;(fi, . . - ft,,) and the convention Cq = 1. We adopt the definition and notations for W and 2V as given 
by (|3.2| | and | |3.3|I respectively, along with this renaming of the independent variables. Let us define the set 
T :- {ti, . . . ,1^} and the set omitting the variable tj by Tj := T\{tj]. Furthermore we define the Vandermonde 
determinant 

(4.35) A(r)- Yl ih-tj). 

l<j<k<N 



The defining relation for the weight l|2.1|l implies that the Toeplitz elements {Wk]'^^_^ satisfy the N+ 1-order 
linear difference equation 



N+l 



(4.36) Yj'^-i [(;■ - /)eN+i-/ - m^+i-i] iVj-i = 0, je Z, 



with N + l consecutive elements being arbitrary "initial" values. These "initial" values also define the Li 
polynomial, which will in turn fix the initial values of our recurrence relations I4.3[ by the expression 



N+l 



(4.37) ^-Yj 



uiz', 



1=0 
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where the coefficients are 

(4.38) Mo = {-)^womN+i, "n+i = woMq, 

i-i 

(4.39) Uj = {-f-'woniN^i.j + 2Y^(-f^'-' [(; - 1)6^^^., - mw+i.,] Wj.,, j=l,...,N. 

1=0 

We adopt a parameterisation of the spectral coefficients of the form 

N-l 

(4.40) — 0„(z) = i-fine^^i - mw+i)— + Y 9'„z' + {n + l+ mo)z^ 
and 

(4.41) Q„(z) = (-)^(«eN+i - -mw+i) + (-)~ ^^^(^ + -(-)'m;v+i_Oz' + (1 + -m,)z''^\ 

1=1 

which introduces the variables d],, . . . , ^JJ'"^ and aj],, . . . , co^. 
Proposition 4.3. Define the variables 

The following system oflN coupled first order recurrence relations in nfor the variables [fl, w^j^j 
[{n - po) Ylk^j h + Lti t'r^c^'n + {-fa + mo)tf] 



(4.43) tjflfi^^ 



[in - po) Uk ik + Lti 4 + (-)~(1 + mo)] 



[n Uk*j ik + Llr t'r\a,'„ + (-)'«%,!_,) + (-)~ff 

[n Uk h + LlMr + i-yniN^i-i) + i-r] 



and 

(4.44) co'„ + col_^ + (-)%N+i-y = {-y{n - l)eN+i-;(r U {!)) 



A{T)eN-j{T) + Lti{-f-'%A{TiU{l})eN.j{TiU{l])f!, 
+ {-Yin - po) 



A(r) + i:f^i(-)~-i-'A(r,u{ii)/i 

A(r)e^„i_/r) + i:r=i(-)~-^"'f/A(T/U{l|)e^,,i_/T,U {!))/,( , 

(-)^(n + l + mo) ■ T, ■ / / = 1,---,N, 

A(T) + Lf.i(-)~-i-'f,A(r;U|l})/i 



completely characterises the bi-orthogonal polynomial system. The recurrence relations are subject to the initial 
conditions 



(4.45) fi 



_ 2V(tj)-KlU (tj) 
f .m//1^ _ 1.^27 



7 _ 



tj[2V{l) - kIU{1)] 

(4.46) (-)^a;^ = -^[z'] (2V + k^LT) - ^[z^-^] (2V - k^LI) , 

/or = 1, . . . , N, ivhere [z'](.) denotes the coefficient of J in the polynomial. 
Proof. To establish (|4.431 we utilise l|2.18t in the following form 

tj@n{tj)@n^l{tj) ^ [nSj) + V^(i;)][n«(t;) - Vitj)] 

^' ' e„(i)e„+i(i) [Q„(i) + y(i)][Q„(i)-y(i)] ' ^ 
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We note that the factors appearing on the right-hand side can be found from 



(4.48) 



(4.49) 



1=1 

N 

a -V= {-fne^^i + (-f Yj{J„ + (-)'mN+i-Oz' + z^^^ 



Then l|4.431 follows from this result and the definition ||4.42| |. 

To prove the relations ||4.44| | we first need to invert the definition | |4.42| | along with the parameterisation 
(|4.40l l for the coefficients r„/r„+i and di, f or ;■ = 1, . . . , N. We find 

, . _ , , A(T) + Lt^{-r-'''A{T, u ml 

(4.50) (n-po) = (n + l+mQ)- 



rn+i 



and 
(4.51) 



9'n = {-f^^n + 1 + mo) 



A{T)eN-j{T) + Lli{-f-'%HTi U {l])eN-i{T, U {!])/'„ 



A(T) + Lli{-r-'^'t,A(T, U {1})/^ 
for = 1, . . . , N - 1 by using the identity for Vandermonde determinants Aj(T) with the ;'-th column missing 



(4.52) 



Aj{T) := det 



(1 h 



1 



[] 



[] 



eN-;(T)AN(T), 



for i = 0,. . .,N. The initial values of the dependent variables (|4.45l and ||4.46|I are found using (|2.231 and 
(|2.25t respectively and their definitions. □ 

Remark 4.1. Of primary interest in applications are the Toeplitz determinants 1 11. 21 which are also T-fimctions 
in the sense of Jimbo-Miwa-Ueno's definition Il23ll , l2l1 ,||22 | . Recovery of the x-function /„ from the f,{, o)/, of 
the previous proposition can be achieved by solving the recurrence in r„ using | |4.50|I and then by employing 
either of 



(4.53) 
or 

(4.54) 



1 = -(n + l)ei - mi + (m + 2 + mo) (— - A„+2 ) + {n + mo)A„, 
i-fa)^ = -ei-mi + (n + l+ mo)A„+i + (n + 2 + mo) (— - An+^, 



and solving these recurrences for A„ . This latter sequence can then used to find r„ via l|1.9l l, and with both 
of the r-coefficients one can use l ll.lOI I to solve for the recurrence in /„. 

Remark 4.2. The mapping between the Hamiltonian variables for and those of the above proposition are 
given by l|4.42ll with (|3l4t and 



(4.55) 



- mqr)Vr = {-f{n - po)eN.i + i-lfY,coU + (1 + mo)qf^\ 



or 



(4.56) 



1 + mo + (n - po) 



CNiT) 



eN{Q) 



En- 



-;(Q)-£eN-y(a) 



)-l 



icir-l)nl,{qr-k) 



where Q := l^sl^j and Qr := QWqA- Written in these Hamiltonian co-ordinates the second of the coupled 
recurrences | |4.44|| is 



(4.57) 



_ n 2V{qr) 

Vn+l,r + Pn,r — 777} — T ■ 
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Remark 4.3. A system of discrete Gamier equations could be formulated using the "conjugate" variables 0* 
and Q* but this would not differ in essence from the one presented here. 

Remark 4.4. We have not attempted to check whether singularity confinement |12|, an algebraic entropy 
criteria IITSi , ID, ||35| or one based on Nevalinna theory ||2l, Il4i applies to our recurrences and this remains 
an outstanding issue. 

Remark 4.5. There have been reports of recurrence relations for the Gamier systems in fSSl and fSTl, however 
the explicit relationship between these equations and the recurrences reported here remains to be elucidated. 

Remark 4.6. Another possible method for deriving the recurrence relations of Proposition 14.31 would be to 
construct the Schlesinger transformation, or lattice translation, operators from the fundamental reflection 

and automorphism operators of the affine Weyl group B^Jl^g, which have been studied in 1 26 1, 133, 1331 . This 
task was carried out for M = 3,N = 1, i.e. for the sixth Painleve equation in |8|, however this involved a 
laborious calculation and it may not be feasible to employ this approach to the multi-variable extension. 

In this study we have focused on a particular type of transformation, namely that of n i-^ n ± 1 or 
equivalently do t-^ do ± 1 and 0oo ^ 6oo ± 1, which is natural within this context. However one could 
legitimately ask for the recurrence systems for the transformations py i-> py ± 1, which are part of the larger 
group of Schlesinger transformations. The theory of these, in the context of bi-orthogonal system on the 
unit circle, has been investigated in [351 but analogues of the recurrences found here were not given there. 
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